Absence of a structural glass phase in a monoatomic model liquid predicted to 

undergo an ideal glass transition 



o 
o 

C/3 



I 



£ 

i 

C 

O 
o 



> 

m 
m 

co 

o 
o 



Charlotte Gils, Helmut G. Katzgraber, and Matthias Troyer 
Theoretische Physik, ETH Zurich, CH-8093 Zurich, Switzerland 
(Dated: February 1, 2008) 

We study numerically a two-dimensional monodisperse model of interacting classical particles 
predicted to exhibit a static liquid-glass transition. Using a dynamical Monte Carlo method we 
show that the model does not freeze into a glassy phase at low temperatures. Instead, depending 
on the choice of the hard-core radius for the particles, the system either collapses trivially or a 
polycrystalline hexagonal structure emerges. 

PACS numbers: 64.70.Pf, 61.43. Fs, 61.43.Bn 



I. INTRODUCTION 

There has been considerable progress in the physics of 
structural glasses in the last few years. 1-2 While experi- 
mentally as well as numerically new insights have been 
gained, theoretical studies of simple glass-forming mod- 
els lag behind. A structural glass is a solid state char- 
acterized by 'random' particle positions such that it is 
not possible to identify symmetries, i.e., no long-range 
order exists. The nature of the transition from a liquid 
to such a disordered phase is still an open question (see 
for example Refs. and0). A popular microscopic ap- 
proach is that of mode-coupling theory which predicts a 
purely dynamical transition, i.e., a change in the dynam- 
ics from ergodic to non-ergodic behavior (see for example 
Ref. [H). In contrast, a static glass transition arises from 
an equilibrium classical statistical mechanics approach 
using a replica formulation^^ based on the assumption 
of an entropy crisis scenario.? It is important to note 
that this approach assumes that the system is in a uni- 
form phase, thereby ignoring a crystal state from the very 
outset. A static glass transition may therefore only arise 
in systems where a crystal state does not dominate the 
low-temperature phase, e.g., by being separated from the 
glass phase by a large energy barrier. 

The replica technique has been applied to a variety of 
systems, such as identical hard spheres^ a Lennard-Jones 
binary mixture^ and to a monoatomic model with at- 
tractive two-body interactions*^ The last study has the 
advantage over conventional glass-forming models in that 
exact replica calculations can be performed thus poten- 
tially being able to deliver new insights into this com- 
plex field. The authors of Ref. 0j] consider a system of 
N classical particles in D dimensions which interact by a 
two-body potential which is nonzero for a given distance 
R between the particles and zero otherwise (see Fig. [T|). 

Analytica l 12 : 13 and numerical studies^ report the ab- 
sence of a static glass phase in monodisperse hard spheres 
in two and three dimensions, respectively— Therefore 
it would be of interest to test if the model proposed 
in Ref. [ll| has a static glass transition even though the 
model is intrinsically monodisperse [see Eq. ([2|) below] . In 
this work we study the model proposed by Dotsenko and 
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1: Schematic plot of the interaction potential U(r) 
©] between two particles with relative distance r = 
For r £ [R — ro,R + ro] the particles feel an attrac- 
tion, otherwise the potential is zero. Note that the model is 
defined in the limit 1 <C ro <C R. (a) Interaction potential as 
introduced in Ref. |Tl] with a hardcore radius h = 1. (b) Mod- 
ified interaction potential in order to avoid a trivial particle 
collapse with a hardcore radius R/4 < h < 72/2. 



Blatter— numerically in two space dimensions (D = 2) 
and find that for various combinations of the model pa- 
rameters the system crystallizes immediately. We thus 
conclude that some amount of disorder must be intro- 
duced in order to stabilize the glass phase. 



II. MODEL 

We study the monoatomic glass former proposed in 
Ref. [ll]. The model describes N particles in D space 
dimensions with a two-body interaction potential 
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for i?_ < \t{ — Tj\ < R+. Here R± = R ± ro, where R is 
the 'interaction radius' of the potential and ro the width 
of the parabolic potential well [see Fig.[T]for details]. The 
particle size is set to unity, which is equivalent to a hard- 
core radius of h — 1. In these units the results obtained 
in a mean-field approximation are valid in the parameter 
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FIG. 2: Configuration of particles after 2 • 10 8 Monte Carlo 
steps for r = 10, R = 100, h = 30 (i.e., T c w 0.33), T = 0.1, 
N = 120 and an average particle spacing p _1//2 ~ 88.3. 

regimes 1 <C r <C -R and for an average particle distance 
which equals approximately those of the crystal and the 
liquid phase, i.e. p~ x l D ss R (p denotes the density). The 
liquid and glass phases are predicted to be separated"! 
by a critical temperature 

T c ~l/ln(Die/r ). (3) 

The disordered nature of the phase for temperatures be- 
low T c is identified by means of a replica correlator which 
quantifies the correlations between particles in two or 
more replicas. 

III. NUMERICAL METHOD AND RESULTS 

We use a Metropolis Monte Carlo (MC) method (see 
for example Ref. Ila for details). Without loss of gen- 
erality we study particles in two space dimensions with 
periodic boundary conditions in a box of size L x x L y . A 
MC step consists of choosing a new configuration of par- 
ticle positions and accepting this move with probability 
p = min[l, exp(— AU/T)], where AU is the difference in 
the potential energy between the old and new configura- 
tion of particles. Two types of local moves are applied: a 
random particle and a random direction are chosen, and 
the particle is moved into this direction by a distance 
uniformly drawn from either the interval [0,ro], or the 
interval [R_,R + ]. The resulting local dynamics serves 
as a model for the realistic dynamics of the system in 
renormalized physical time. In the absence of ergodicity 
breaking, the method also allows one to study the ther- 
modynamic phase of the system: after a sufficient num- 
ber of MC steps configurations are distributed accord- 
ing to the equilibrium distribution exp(— U/T). Clearly, 
if no disordered phase is found after relatively few MC 
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FIG. 3: Configuration snapshot after 3 • 10 9 Monte Carlo 
steps for ro = 10, R = 100, h = 30, (i.e., T c « 0.33), T = 0.1, 
N = 1200 and an average particle distance p _1//2 ~ 108 (i.e., 
a smaller density than in Fig. [2}. A polycrystalline structure 
with randomly-oriented patches of hexagonal crystal is found. 



steps, the use of optimized nonlocal MC updates which 
allow the system to assume low-energy configurations in 
a much faster manner will certainly not lead to a disor- 
dered phase either. The simulation is started with a ran- 
dom configuration of particles. The temperature T for 
the simulation is chosen such that T < T c , where T c is 
given by Eq. ([3|) . Since the system starts from a quenched 
configuration, the cooling rate is thus infinitely high, as 
required in Ref. [Til 

For the glass phase to be stable Dotscnko and Blatter— 
give clear restrictions on the parameters of the model 
glass former, i.e., 1 <C ro <C R. Simulating the system 
with the standard hard-core radius h = 1 and for temper- 
atures T <T C yields a trivial collapse where all particles 
approach each other up to a minimum distance of 2h. We 
have tested this scenario for various choices of all system 
parameters. The same trivial collapse is found for any 
hard-core radius which satisfies h < R/4. To avoid this 
behavior, we increase the size of the hardcore radius to 
values R/2 > h > R/4. Note that the size of the hard- 
code radius h is irrelevant in the analytic calculations of 
Ref.HTl.i2 

For the modified interaction potential [Fig. [IJb)] a 
hexagonal crystal structure emerges: If L x , L y and N are 
chosen appropriately, a perfect crystal develops as shown 
in Fig. [2j If the density is lower than the density of the 
hexagonal crystal, the particles organize in a polycrys- 
talline structure: small hexagonal crystal patches with 
random orientations form, as can be seen in Fig. [3] and 
Fig. Hid). In Figs. Ufa) ~2ld) we show how the crystalline 
structure forms as a function of time measured in Monte 
Carlo steps, starting from a low-density random config- 
uration [panel (a)]. We have also measured the replica 
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FIG. 4: Particle configurations for N = 1000 particles (r = 10, R = 100, h = 30, T = 0.05 < T c ) in a system with 
average particle distance p _1//2 ~ 167. (a) Initial random configuration, (b) configuration snapshot after 2 ■ 10 7 MC steps, (c) 
configuration snapshot after 2 • 10 8 MC steps, and (d) after 5 ■ 10 9 MC steps. As in Fig. [3] a polycrystalline structure emerges. 



correlator [Eq. (62) in Ref. ED] whose finite value would 
be indicative of a glassy state. However, in our simula- 
tions, it takes a large value close to the linear system size 
squared, which means that particles in two replicas are 
completely uncorrelated (not shown). For large systems 
and various densities, ranging from those of the crystal 
state [Fig. [5] to much more dilute systems [Figs. [3] and 
[4], our results arc in stark contrast to the predictions 
of Ref. [ll|: While a glassy "network-like" structure with 
an average number of nearest-neighbors which is signifi- 
cantly smaller than in a crystalline state is predicted, our 
results are locally crystal-like and show a high degree of 
orde n 18 i 19 We have varied the system parameters, in par- 
ticular increasing the ratio of R and ro, but the results 
have been qualitatively similar to the ones shown in the 
Figures. Therefore, if a glass phase exists in the model, it 
cannot be observed since crystallization takes place very 
quickly. 



body interactions, such as the recent proposal^ studied 
here in detail, cannot form a glassy structure at low tem- 
peratures. Instead, the system forms a polycrystalline 
structure and does not show an ideal liquid-glass transi- 
tion. This is unfortunate since the proposed monoatomic 
model is analytically solvable and thus holds the promise 
to deliver new insights to the physics of glass formers. 
It might be of interest to verify if the model is still ana- 
lytically solvable if, for example, many-body interactions 
between the particles 2 ^ are added to the Hamiltonian or 
Gaussian randomness is introduced into the distance R 
to suppress the crystal state and favor a glassy phase. 
Assuming the crystal structure is better suppressed in 
higher space dimensions, it might be conceivable to find 
a glass structure for D > 2 and thus possibly determine 
a lower critical dimension of the model. 



IV. CONCLUSIONS 

In conclusion, our results from Monte Carlo simula- 
tions are in agreement with previous prediction a 12 ' 13 that 
monoatomic systems in two space dimensions with two- 



While the simulations are performed for finite systems, 
we believe the number of particles studied shall not in- 
fluence the conclusions found. In particular, we have 
scanned the vast parameter space in detail and have per- 
formed simulations for many different particle numbers 
and densities always yielding similar results. 
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